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During its first six years of operation, the Fermi Large Area Telescope (LAT) has detected >30 MeV gamma- 
ray emission from more than 40 solar flares, nearly a factor of 10 more than those detected by EGRET. These 
include detections of impulsive and sustained emissions, extending up to ~20 hours in the case of the 2012 
March 7 X-class flares. We will present an overview of solar flare detections with LAT, highlighting recent 
results and surprising features, including the detection of >100 MeV emission associated with flares located 
behind the limb. Such flares may shed new light on the relationship between the sites of particle acceleration 
and gamma-ray emission. 


1. Introduction 

Understanding the processes of particle accelera¬ 
tion and impulsive energy release which occur in nu¬ 
merous sites throughout the Universe is one of the 
major goals of space physics and astrophysics. The 
Sun is the most powerful particle accelerator in the 
solar system and its proximity permits investigat¬ 
ing the entire electromagnetic spectrum of these ac¬ 
celeration phenomena. During solar flares, the Sun 
is capable of accelerating electrons and ions to rel¬ 
ativistic energies on time scales as short as a few 
seconds, as indicated by observations of X-rays, mi¬ 
crowaves, 7-rays, and neutrons produced when the 
flare-accel erated particles interact with the solar at - 
mosphere [Forrest and ChuppIllQS.'lllKane et al.lll98dj |. 
In general, the 7-ray emission light curve is similar to 
that of the HXRs (possibly with some delay), last¬ 
ing for 10-100 seconds. This is referred to as the 
“impulsive” phase of the flare. However, the En¬ 
ergetic Gamma Ray Experiment Telescope (EGRET) 
on-board CGRO [Kanbach et al.lilQ^ lEsposito et al.l 
Il999l| also detected a sustained emission in gamma 
rays for m ore than an h our after the impulsive phases 
of 3 flares [Rvanll^OOaj . The expected increase of so¬ 
lar activity during the current solar maximum is pro¬ 
ducing a large number of solar flares, including bright 
GOES X-class and moderate M-class flares. 


sitive from ~ 8 keV up to 40 MeV [Meeean et al.l 
|2009|| . During the first 18 months of operation co¬ 
inciding with the solar cycle minimum, the Fermi 
LAT detected >10 0 MeV gamma-ray emission from 
the quiescent Sun [Abdo etldl l2011a |. As the so¬ 
lar cycle approaches it maximum, the LAT has de¬ 
tected several solar flares above 30 MeV during both 
the impulsive and the temporally ex t ended phases 
JOhnoet_aL 2011. Omodei et al. 2011. ^naka et ^ 


I2OI2L ^trosian et al.l 201^ Omodei et al.l2012| . The 

first Fermi GBM and LAT detection of the impul¬ 
sive GOES M2.0 flare of 2010 June 12 is presented 
in lAckermann et al] [2012a |. The analysis of this 
flare was p erformed using the LAT Low-Energy (LLE) 
technique Pelassa et al.ll201(]j| because the soft X-rays 
emitted during the prompt emission of a flare pene¬ 
trate the anti-coincidence detector (ACD) of the LAT 
causing a pile-up effect which can result in a signifi¬ 
cant decrease in gamma-ray detection e fficiency in the 
standa rd on-ground photon analysis [Atwood et al.l 
l2009al . The pile-up e f fect ha s bee n addressed i n detai l 


Ackermann et al. I [2012a| and lAbdo et al.l [20n9l| . 


The list of all LAT detected flares, and the analysis of 
the first two flares with long lasting high-energy emis- 
si on (2011 March 7- 8 and 2011 June 7) is pr esented 
in lAckermann et al] [2014 1 , lAiello et al.l [2014 j. 


2. Fermi observations of the Sun 


2.1. June 2010: An impuisive event 


Launched in 2008, the Fermi observatory is com¬ 
prised of two instruments; the Large Area Telescope 
(LAT) designed to detect gamma-rays from 20 MeV 
up to more than 300 MeV Atwood et al.ll200'^ and 
the Gamma-ray Burst Monitor (GBM) which is sen- 


On 2010 June 12 00:30 UT a moderate GOES M2.Q 
class X-ray flare erupted from the active region (AR) 
11081 located approximately N23°W43°. At the time 
of the flare the Fermi spacecraft was in sunlight and 
during a relatively low-background portion of its or- 
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bit^. The GBM triggered on the flare at 00:55:05.64 
UT and detected keV emission for RilO m. The 11 — 26 
keV emission recorded by the GBM Nal detectors rose 
precipitously for about 40 s and is shown in Figure[T^. 
For comparison we include the GOES 0.5 - 4 Aprofile 
and note that this emission is dominated by 3 keV 
thermal photons as is reflected in its slower rise and 
extended tail. The 100 — 300 keV time profile observed 
by the GBM’s solar facing Nal detector is also plotted 
in Figure [1^. It is clear that the emission peaks more 
sharply and ends sooner at higher X-ray energies. 

The accompanying hard X-ray emission from the 
flare was detected in the LAT’s ACD and is reflected in 
the shape of the average number of AGD tile hits as a 
function of time (shown in Figure[T})). The broad peak 
with a maximum near 00:57 UT of the hit distribution 
has a shape similar to the 11-26 keV emission and 
the impulsive peak is similar to the 100 to 300 keV 
flux observed by the GBM Nal detector. As shown by 
the red curve in Figure [TJ: there is no evidence for the 
flare in the well-screened standard LAT data products. 
[What is shown here are the events belonging to the 
P6TRANSIENT event class. lAtwood etliribOOQBl . This 
is the direct consequence of the pulse pile-up effect. 
The black curve in Figure (TJ; is the LAT LLE >30 
MeV event rate for the time of the flare. 


White light emission observed by the Helioseismic 
and Magnetic Imag er (HMI) on the Solar Dynam ics 
Observatory {SDO) [Martinez Oliveros et al.lf20II in 
a single 45 s exposure during the hard X-ray emission 
revealed two compact footpoints about 10'^ km apart. 


The >30 MeV LLE spectrum of this flare revealed 
flare emission up to an energy of ^400 MeV. The 
nuclear line emission observed with the GBM im¬ 
plies the presence of accelerated ions up to at least 
50 MeV nucleon”^. It is possible that the flare- 
accelerated proton spectrum extended up to the ~300 
MeV threshold for pion production. Alternatively, it 
is also possible that the LAT emission is from elec¬ 
tron bremsstrahlung, either from an extension to high 
energies of the electron spectrum producing the X-ray 
bremsstrahlung observed in the GBM or from an addi¬ 
tional hard electron component. One possible way to 
resolve this ambiguity is to jointly fit the GBM and 
LAT spectra assuming different origins for the LAT 
emission. 


In Figure[2]we plot the background-subtracted pho¬ 
ton spectrum from 0.3 to 400 MeV including both the 
GBM and LAT data. We made two fits, using rmfit 
3-4^, customized for the specific solar flare, and the 


^The Fermi observatory is in a nearly circular orbit with an 
inclination of 25.6° at 565 km. 

^R.S. Mallozzi, R.D. Preece, & M.S. Briggs, “RM¬ 
FIT, A Lightcurve and Spectral Analysis Tool”, Robert 
D. Preece, University of Alabama in Huntsville, (2008): 


OSPEX^ analysis packages, to the joint data sets. In 
the first fit we assume that the observed LAT emis¬ 
sion was from pion-decay radiation (top panel of Fig¬ 
ure and the other assuming that it was from a hard 
power-law spectrum of electron bremsstrahlung (bot¬ 
tom panel). Based on the statistical quality of the 
fits to the LAT spectrum we cannot distinguish be¬ 
tween the two emission models but, if the LAT emis¬ 
sion is from electron bremsstrahlung, we have found 
that it cannot be a simple extension of the low-energy 
bremsstrahlung components that we determined from 
fits to the GBM data; it must be from a distinct popu¬ 
lation of electrons extending to energies of several hun¬ 
dred MeV. However, this high energy electron compo¬ 
nent would produce a spectrum that steepens beyond 
tens of MeV due to syn chrotron energy losses that 
increase with energy [see lPark et al.l[l997| |. and must 
have a quite different origin. Gonsequently we believe 
that this is a less likely scenario than the hadronic 
model. 


Assuming that the LAT emission is from hadronic 
interactions, we have fit the L AT spectrum with cal - 
culated pion-decay templates [Murphy et al.l Il987 |. 
which depends on the ambient density, composition 
and magnetic field, on the accelerated-particle compo¬ 
sition, pitch angle distribution and energy spectrum. 
The templates represent a particle population with an 
isotropic pitch angle distribution and a power-law en¬ 
ergy spectrum {dN/dE oc with E the kinetic en¬ 
ergy of the protons) int eracting i n a th ick target with 
a coronal composition [Reame d Il995l| taking ^He/H 
= 0.1. With 67% confidence (based on x^) we con¬ 
clude that the spectrum of accelerated ions responsi¬ 
ble for the pion-decay emission must be steeper than 
a power-law with index —4.5. We note that there is 
no change in the quality of the fits for indices steeper 
than —5 due to limited statistics >400 MeV. We can 
use the results of our GBM and LAT spectral anal¬ 
yses to obtain information on ions acc elerated in the 
impul sive phase of the June 12 flare. iMurphv et al.l 
[l997 | have described how parameters derived from 
integrated spectroscopic fits and temporal studies can 
be used to obtain this information. We first use the 
nuclear de-excitation line, neutron-capture line, and 
pion-decay fluences to estimate the overall shape of 
the accelerated ion spectrum. These three emissions 
are produced by accelerated ions within distinct en¬ 
ergy ranges: ~5-20 MeV for the de-excitation lines, 
~10-50 MeV for the neutron capture line, and >300 
MeV for the pion-decay emission. Ratios of these 
emissions therefore determine the relative numbers of 
accelerated ions in the associated energy ranges. We 
then obtain spectral indices across these energy ranges 


http: / / fermi.gsfc.nasa.gov/ssc / data/analysis / user/ 
'“SolarSoft: http://www.lmsal.com/solarsoft/ 
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Figure 1: Time histories related to the 2010 June 12 solar flare, a) GOES 0.5 - 4 Arates, and GBM Nal 11 - 26 keV 
and 100 - 300 keV relative rates; b) LAT ACD hit rate >100 keV containing contributions from background, >100 keV 
solar flare X rays (impulsive peak) and pulse pile up from lO’s of keV solar X rays following the Nal 11 - 26 keV profile 
in la); and c) LLE and LAT Transient Class event rates. 


by comparing meas ured ratios with ratios from theo- 
retical calculations Murohv et al.l Il987l l2005l l2007j 
based on updated nuclear cross sections. 

If we assume that the LAT emission >30 MeV was 
entirely due to pion-decay emission, then we estimate 
that the flare-accelerated ion spectrum was consis¬ 
tent with a series of power laws, softening with en¬ 
ergy, with indices of ^—3.2 between ~ 5 — 50 MeV, 
4.3 between ^50-300 MeV, and softer than ^—4.5 
above 300 MeV. In Table U we summarize our find¬ 
ings, reporting the processes responsible for the de¬ 
tected emission, energy range of emitted gamma-rays, 
as well as the energy and spectral index of the accel¬ 
erated ions/electron distribution. 


Component 

7 -rays 

electrons/ions Spectral Index 


(MeV) 

(MeV) 

acc. particles 

Brem. 

0.1-1 

0.1-1 

-3.2 

Brem. 

2-10 

2-10 

<-1.2 

HE Brem. 

10-200 

10-200 

Ri-2.0 

Neutron Capt. 

2.2 

5-50 

-3.2 

Nuclear lines 

5-20 

50-300 

-4.3 

Pions 

>30 

A. 300 

- -4.5 


Table I Derived quantities for accelerated particle 
distributions (with a cut-off at 2.4 MeV) 


2.2. March 2012: Impulsive and 
sustained emission of a bright fiare 

On 2012 March 7 two bright X-class flares origi¬ 
nating from the AR NOAA AR^: 11429 (located at 
N16°E30°) erupted within an hour of each other, 
marking one of the most active days of Solar Cycle 


24. The first flare started at 00:02:00 UT and reached 
its maximum intensity (X5.4) at 00:24:00 UT while 
the second XI.3 class flare occurred at 01:05:00 UT, 
reaching its maximum 9 minutes later. 

The GOES satellite observed intense X-ray emis¬ 
sion beginning at about 00:05:00 UT and lasting for 
several hours. Moreover, it detected Solar Energetic 
Particles (SEP) protons in three energy bands origi- 
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Figure 2: Combined GBM/LAT photon spectrum 
accumulated between 00:55:40 and 00:56:30 showing the 
best total fit using the same components as in Figure 3 
plus an additional component for the LAT emission. The 
upper panel shows a pion-decay fit to the LAT spectrum; 
alternatively the lower panel shows a power-law fit, 
presumedly representing a third electron bremsstrahlung 
component. Note that because this is a photon 
representation the lines are plotted at their intrinsic 
resolution and appear to be more significant than they 
really are. 


nating these flares. The Reuven Ram aty High-Energy 
Solar Spectroscopic Imager [RHESSI, iLin et alll200^ 
was not obserying the Sun during this period. On the 
top panel of Figure |3] we plot the X-ray data from 
GOES 15 satellite measured in both 3-25 keV and 
1.5-12 keV, as well as the detected proton flux. 

The Fermi LAT >100 MeV count rate was dom¬ 
inated by the gamma-ray emission from the Sun^, 
which was nearly 100 times brighter than the Vela 
Pulsar in the same energy range. During the impul¬ 
sive phase (the first eighty minutes) the X5.4 flare was 
so intense that the LAT observation suffered from the 
pile-up effect so we used the LLE technique to analyze 
the impulsive phase of this bright flare. 


^http: //apod.nasa.gov/apod/apl20315.html 


We fit the data using XSPEC^ to test three models. 
The first two are simple phenomenological functions, 
to describe bremsstrahlung emission from accelerated 
electrons, namely a pure power law (PL) and a power 
law with an exponential cut-off (EXP): 


dN(E) 

dE 


= Nqc ^ exp 



( 1 ) 


where P is the photon index and E^o is the cut¬ 
off energy. We found that the data clearly diverge 
from a pure power law spectrum and that the EXP 
provides a better fit in all time intervals considered. 
The third model use d the same pion decay templates 
Murphy et ahlllQ^ I used for the 2010 June 2 flare. 


When using the pion-decay templates to obtain the 
gamma-ray flux value we fit the data varying the pro¬ 
ton spectral index from 2 to 6, in steps of 0.1. In this 
way, we fit the LAT data with a model with two free 
parameters, the normalization and the proton index 
s. 


To study the temporally-extended emission, we per¬ 
form time-resolved spectral analysis in Sun-centered 
coordinates by transforming the reference system from 
celestial coordinates to ecliptic Sun-centered coordi¬ 
nates. This is necessary in order to compensate for 
the effect of the apparent motion of the Sun during 
the long duration of the flare. We select intervals when 
the Sun was in the FOV (angular distance from the 
LAT boresight < 70°) and use the unbinned maximum 
likelihood algorithm implemented in gtlike^. 

We include the isotropic template model that is 
used to describe the extragalactic gamma-ray emis¬ 
sion and the residual cosmic ray (CR) contamina¬ 
tion®, leaving its normalization as the free parame¬ 
ter. Over short time scales, the diffuse Galactic emis¬ 
sions produced by GR interacting with the interstellar 
medium are not spatially resolved and are hence in¬ 
cluded in the isotropic template. We also add the 
gamma-ray emission from the quiescent Sun modeled 
as a point source located at the center of the disk, 
with a spectrum described by a simple power law with 
a spectral index of 2.11 and an integrated energy flux 
(> 100 MeV) of 4.7x10“^° ergcm“^s“^ [correspond- 
in g to a flux of 4.6x l0~^ phcm“^s“^ as reported 
in lAbdo et al.l l201l'9| . We did not include the ex- 
te nded Inverse Comp ton (IC) component described 
in lAbdo et al.l [2011b| because it is too faint to be de¬ 
tected during these time intervals. We fit the data 
with the same two phenomenological functions used 


^http: //heasarc.gsfc.nasa.gov/docs/xanadu/xspec/index.html 
“^We used ScienceTools version 09-28-00 avail¬ 
able from the Fermi Science Support Center 
http://fermi.gsfc.nasa.gov/ssc/ 

^We used iso_p7v6source.txt available from the Fermi Sci¬ 
ence Support Center 
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Figure 3: Long lasting emission. Top panel: soft X-rays (red: 1.5-12 keV, blue: 3-25 keV) from the GOES 15 satellite. 
On the right axis, 5-minute averaged proton flux (green: 30-50 MeV, yellow: 50-100 MeV, magenta: >100 MeV). We 
display the average of detectors A and B. Bottom panel: high-energy gamma-ray flux above 100 MeV measured by 
the Fermi LAT. The Blue/red circles represent the flux and the derived proton spectral index obtained with the LLE 
analysis (covering the initial period, when the instrumental performance was affected by pileup of hard-X-rays in the 
ACD tiles). The blue circles and red squares represent the flux and derived proton spectral index, respectively, 
obtained by standard likelihood analysis. Green diamonds are the GOES proton spectral indexes derived from the 
hardness ratio, as described in the text. The gray bands correspond to the systematic uncertainty associated with flux 
measurements and of the estimated proton index due to uncertainties on the effective area of the i nstrument. The 
horizontal dashed line corresponds to the value of the gamma-ray flux from the quiescent Sun, from lAbdo et ah] |2011b|| . 


for the impulsive phase of the flare and use the like¬ 
lihood ratio test to estimate whether the addition of 
the exponential cut -off is statistically sig nificant. The 
Test Statistic (TS) iMattox et al.l (l996l | is defined as 
twice the increment of the logarithm of the likelihood 
C obtained by fitting the data with the source and 
background model components simultaneously. Be¬ 
cause the null hypothesis (i.e. the model without an 
additional source) is the same for the two models, the 
increment of the TS (ATS=TSplexp-TSpl) is equiv¬ 
alent to the corresponding difference of maximum like¬ 
lihoods computed between the two models. 

For each interval, if ATS > 30 (roughly correspond¬ 
ing to 5cr) then the PLEXP model provides a signif¬ 
icantly better fit than the simple power-law and we 
retain the additional spectral component. In these 
time intervals, we also used the pion decay model to 
fit the data and estimated the corresponding proton 


spectral index. We performed a series of fits with the 
pion decay template models calculated for a range of 
proton spectral indices. We then fit the resulting pro¬ 
file of the log-likelihood function with a parabola and 
determine its minimum (£min) and the correspond¬ 
ing value So as the maximum likelihood value of the 
proton index. 

In the lower panel of Figure [3] we combine the LLE 
and likelihood analysis results, showing the evolution 
of both the gamma-ray flux and the derived spectral 
index of the protons®. In the last five time intervals 
the power-law representation is adequate to describe 


®After approximately 11:00:00 UTC the flux of the Sun di¬ 
minished to the point that the spectral index of the proton 
distribution cannot be significantly constrained. 
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the data; in the last bin, the flare is only marginally 
significant (TS=7); the flux and the photon index are 
compatible with the values of the quiescent Sun. For 
this reason we have indicated the last point as an up¬ 
per limit (computing the 95% C.L.). Unlike during 
the impulsive phase, the spectrum during the tempo¬ 
rally extended phase becomes softer monotonically (s 
increases). 

We also compare our results with the GOES proton 
spectral data. For this, we selected two energy bands 
(>30 MeV and >100 MeV) and corrected the light 
curve by the proton time-of-flight (TOF) to 1 AU by 
considering the TOF for 30 MeV and 100 MeV protons 
(i.e. the maximum delay in each energy band). As 
a measure of the spectral index of the SEP protons 
(ssEp)) we compute the Hardness Ratio HRp defined 
as: 

HRp = In ^2) 

P >30MeV 

where P is the integral of the proton flux (assuming 
that the proton flux is proportional to a power law). 
The HRp is related to the value of the spectral index, 
SSEP, of the SEP protons observed at 1 AU, roughly 
as: 


ssEP ^ 1 — 0.83 HRp 


( 3 ) 


To estimate the uncertainty associated with this 
procedure we repeat the calculation neglecting the 
TOF correction. In this way we obtain two values 
for the SEP spectral index for each time bin, corre¬ 
sponding to the actual and zero delay due to the time 
of flight. In Figure [3] we report the estimated pro¬ 
ton spectral index as the average of these two values 
and its uncertainty as half the difference of these two 
values. However we note that the ssep is for protons 
with energy less than a few hundred MeV while s is for 
protons with energies greater than 300 MeV. Diffusion 
is expected to play an important role in the transport 
of these SEPs but an in-depth transport analysis is 
beyond the scope of this paper. From our comparison 
we find that the proton spectral index inferred from 
the gamma-ray data is systematically softer than the 
value of the index derived directly from SEP observa¬ 
tions but that the temporal evolution (hard-to-soft) is 
similar. 

Uncertainties in the calibration of the LAT intro¬ 
duce systematic errors on the measurements. The un¬ 
certainty of the effective area is dominant, and for the 
P7S0URCE_V6 event class it is estimated to be ^10% 
at 100 MeV, decreasing to ^5% at 560 MeV, and in¬ 
creasing to ^10% at 10 GeV and above. We studied 
the effect of the systematic uncertainties on our fi¬ 
nal resul ts via the bracket i ng tech nique described in 
detail in lAckermann et al.l 2012bl| . We find that the 
uncertainties on the flux are < 10% and on the inferred 
proton index are <0.10. The results are represented 
by the gray bands in Figure [S] 


3. New observations: the 
behind-the-limb synopsis 


On 2013 October 11 at 07:01 UT a GOES M1.5 
class flare occurred with soft X-ray emission last¬ 
ing 44 min and peaking at 07:25:00 UT. Figure |4] 
shows the GOES, STEREO-B, RHESSI, Fermi GBM 
and LAT lightcurves of this flare. LAT detected 
>100 MeV emission for ^30 min with a peak flux be¬ 
tween 07:20:00-07:25:00 UT. RHESSI coverage was 
from 07:08:00 — 07:16:40 UT, overlapping with Fermi 
for 9 min. 

Images in Figure [5] from t he STEREO-B Extr eme- 
UltraViolet Imager [EUVI; IWuelser et al.l 12004 1 and 
the SDO Atmos pheric Imaging Assembly [AIA; 
iLemen et al.l l2012j| of the photosphere indicate that 
the AR was ~9.9° behind the limb at the time 
of the flare. LASCO onboard the Solar and He¬ 
liospheric Observatory (SOHO) observed a backside 
asymmetric halo CME associated with this flare be- 
ginnin g at 07:24:10 UT with a linear sp eed of 1200 
km/s [SOHO LASCO GME GATALOGil201^ and a 
bright front over the Northeast. Both STEREO 
spacecrafts detected energetic electrons, protons, and 
heavier ions including helium, as well as type-H 
radio bursts indicating the presence of a coronal- 
heliospheric shock. SDO observed a global EUV wave 
(Liu et al. 2015, in prep.), possibly the coronal coun¬ 
terpart of the shock. STEREO-B had an unblocked 
view of the entire flare and detected^a maximum rate 
of 3.5x10® photons/sec in its 19 5 A channel, cor re- 
sponding to a GOES M4.9 class Nitta et '^l2013|| if 
it had not been occulted. 

The LAT data were analyzed using the unbinned 
maximum likelihood algorithm gtlike implemented 
in the Fermi ScienceTools^ with P7REP_S0URCE_V15 
instrument response functions. We selected gamma- 
rays from a 12° region centered on the Sun and within 
100° of the zenith to reduce contamination from the 
Earth’s limb. For RHESSI data, we implement ed the 
CLEAN imaging algorithm [Hurford et al.l[2002 1 using 
the detectors 3—9 to reconstruct the X-ray images. 
We used the FITS World Coordina te System software 
package [Thompson and \\^l2010l| to co-register the 
flare location between STEREO and SDO images. 
The STEREO light curves are pre-flare background 
subtracted, full-Sun integrated photon rates. 

We measure the direction of the LAT >100 MeV 
gamm a-ray emission [as described in [Aiello et al.l 
[2014 and find a best fit position for the emission 
centroid at heliocentric coordinates of (—855",75") 
with a 68% error radius of 251", as shown in Fig¬ 
ure ISKb). RHESSI X-ray sources integrated over 


^We used version 09-30-01 available from the Fermi Science 
Support Center http://feriiii.gsfc.nasa.gov/ssc/ 
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Figure 4: PRELIMINARY: Light curves of the 2013 
October 11 flare as detected by a) GOES, b) STEREO, 
c) RHESSI, d) GBM, e) LAT and e) RHESSI emission 
centroid heights, with the same color coding as in c). 
Fermi exited the South Atlantic Anomaly (SAA) at 
06:57:00 UT. 


07:11:04—07:16:44 UT are shown as 80%-level, off- 
limb contours in Figure [5l^d). 

The temporal evolution of the projected RHESSI 
source heights above the solar limb are shown in 
Figure lUf). The higher-energy emission generally 
comes from greater heights, cons istent with expecta¬ 
tions for a loop-top source [e.g., iMasuda et alJIlQ^ 
ISui and Holmanl 120031 iLiu et all 2004| . If this were 
a footpoint source, we would expect an opposite 
trend since larger column depths a re required for 
stopping higher-energy el ectrons [e.g.. lLiu et al.ll20^ 
iKontar. E. P. et "^bOOSlI . Moreover, from SDO/AIA 
movies we find no signature of EUV ribbons, even in 
the late phase during the RHESSI night. Together, 
these observations provide convincing evidence that 
the footpoints were indeed occulted. 


4. Discussion 

We have presented the analysis of three solar flares 
detected by the Eermi LAT at high energy, and 



0 200 400 600 800 1000 -1000-800 - 600 -400 -200 0 



500 550 600 650 -950 -900 -850 -800 
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Figure 5: PRELIMINARY: STEREO-B (left) and 
SDO (right) images near the flare peak. The 
white-dashed line in (a) and (c) represents the solar limb 
as seen by SDO. The green line in (b) shows the 68% 
error circle for the LAT source centroid. The cyan 
contour and plus sign in (c) mark the STEREO flare 
ribbon and its centroid, respectively. Their projected 
view as seen from the AIA perspective is shown in (d), in 
which the centroid is located at 9.9° behind-the-limb. 

The green and blue-dotted contours in (d) show RHESSI 
sources. The rectangular brackets in (a) and (b) mark 
the field of view (FOV) for (c) and (d), respectively. 


we highlight some of the similarities and differences 
of these flares. The high-energy emission of the 
2012 June solar flare seems to be correlated with 
HXR emission, suggesting that acceleration of par¬ 
ticles and gamma-ray emission take place close in 
space. Specifically, particles accelerated at the loop 
top could propagate along the loop field lines inter¬ 
acting and emitting gamma-rays at the footprint. For 
this flare, there is no evidence for precipitation of 
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trapped flare particles, particles accelerated in mag¬ 
netic loops after the impulsive phase, particles accel- 
erated in CM E-associated reconnecting current sheets 
[Rvanll2000bj| . or particles sharing the same origin as 
the Solar Energetic Particles (SEPs~) ob served in space 
Ramatv et al.lll98^ ICliver et al.lll993 |. On the other 


hand, flares with long (or sustained) gamma-ray emis¬ 
sion have also been observed by the Fermi LAT. Tem¬ 
poral and spectral analysis suggests that, even if the 
short impulsive phase is clearly visible at > 100 MeV 
energies, the sustained long lasting emission is more 
correlated with SEP properties, suggesting that, for 
this class of flares, either long trapping, continuous 
acceleration, or acceleration at the CME shock could 
be a better explanation. The behind-the-limb hare 
detection at high-energy adds additional considera¬ 
tions that are extremely useful for understanding the 
physics of particle acceleration and gamma-ray pro¬ 
duction during solar flares. We have presented prelim¬ 
inary results from the 2013 October 11 solar flare from 
Fermi, RHESSI, SDO and STEREO. STEREO-R, 
images indicate that the flare occurred in an AR 9.9° 
behind-the-limb. RHESSI and GBM NaR detected 
HXRs up to 50 keV from the flaring loop-top. The 
most unusual aspect of this flare is the LAT detection 
of photons of energies e >100 MeV for about 30 min¬ 
utes with some photons having energies up to several 
GeV. 



Figure 6: Model a): acceleration at the flare, gamma-ray 
emission site below the photosphere; b) acceleration at 
the flare, gamma-ray emission in the corona above the 
limb; c) acceleration (or re-acceleration) at the 
CME-shock, gamma-ray emission at the Sun. 


We consider three scenarios for the emission site of 
the gamma-rays, outlined in Eigure |6l Electrons or 
protons with energies E > e can produce these pho¬ 
tons after traversing a column depth of matter which 
is much larger than the depth penetrated by HXR 
producing electrons (model a). For occulted flares the 
emitted photons must traverse even larger depths than 
particles and they may be scattered and absorbed. Al¬ 


ternatively, acceleration and gamma-ray emission can 
take place in the corona above the limb (model b), 
suggesting trapping of particles, e.g., by strongly con¬ 
verging magnetic fields. In the third model (model 
c) GME-shock accelerated particles can travel back to 
the Sun along magnetic fields connecting the accel¬ 
eration site with the visible side of the Sun. Fermi 
LAT observation of the 2013 October 11 flare (pa¬ 
per in preparation) shows that model a) can be ruled 
out. The LAT detection of gamma-ray emission from 
a flare with 0 > 20° also poses some complications to 
the second scenario (model b), as particles will have 
to be accelerated even further away in the corona, 
where densities are very low. Acceleration (or re¬ 
acceleration) at the CME shock (model c) remains 
possible. Fermi LAT observations are becoming very 
important to disentangle models of particles accelera¬ 
tion and gamma-ray production in solar flares. Future 
LAT observations, combined with a systematic study 
of the solar flares detected at high energy, will very 
likely help to understand this fascinating problem, as 
well as to improve our knowledge of particle accelera¬ 
tion in astrophysical sources in general. 
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